The MetaInvert soil invertebrate genome resource provides insights into below-ground biodiversity and evolution.

Soil invertebrates are among the least understood metazoans on Earth. Thus far, the lack of taxonomically broad and dense genomic resources has made it hard to thoroughly investigate their evolution and ecology. With MetaInvert we provide draft genome assemblies for 232 soil invertebrate species, representing 14 common groups and 94 families. We show that this data substantially extends the taxonomic scope of DNA- or RNA-based taxonomic identification. Moreover, we confirm that theories of genome evolution cannot be generalised across evolutionarily distinct invertebrate groups. The soil invertebrate genomes presented here will support the management of soil biodiversity through molecular monitoring of community composition and function, and the discovery of evolutionary adaptations to the challenges of soil conditions.

Gene abundance linked to climate zone: Parallel evolution of gene content along elevation gradients in lichenized fungi.

Introduction: Intraspecific genomic variability affects a species' adaptive potential toward climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many of these species have remarkable environmental tolerances, and often form populations across different climate zones. Here, we combine comparative and population genomics to assess the presence and absence of genes in high and low elevation genomes of two lichenized fungi of the genus Umbilicaria. Methods: The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked, which genes increase or decrease in abundance within populations along elevation. Results: In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinase superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents and found that the number of this protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations. Discussion: We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.

Metatranscriptomics reveals contrasting effects of elevation on the activity of bacteria and bacterial viruses in soil.

Soil microbial diversity affects ecosystem functioning and global biogeochemical cycles. Soil bacterial communities catalyse a diversity of biogeochemical reactions and have thus sparked considerable scientific interest. One driver of bacterial community dynamics in natural ecosystems has so far been largely neglected: the predator-prey interactions between bacterial viruses (bacteriophages) and bacteria. To generate ground level knowledge on environmental drivers of these particular predator-prey dynamics, we propose an activity-based ecological framework to simultaneous capture community dynamics of bacteria and bacteriophages in soils. An ecological framework and specifically the analyses of community dynamics across latitudinal and elevational gradients have been widely used in ecology to understand community-wide responses of innumerable taxa to environmental change, in particular to climate. Here, we tested the hypothesis that the activity of bacteria and bacteriophages codeclines across an elevational gradient. We used metatranscriptomics to investigate bacterial and bacteriophage activity patterns at five sites across 400 elevational metres in the Swiss Alps in 2015 and 2017. We found that metabolic activity (transcription levels) of bacteria declined significantly with increasing elevation, but activity of bacteriophages did not. We showed that bacteriophages are consistently active in soil along the entire gradient, making bacteriophage activity patterns divergent from that of their putative bacterial prey. Future efforts will be necessary to link the environment-activity relationship to predator-prey dynamics, and to understand the magnitude of viral contributions to carbon, nitrogen and phosphorus cycling when infection causes bacterial cell death, a process that may represent an overlooked component of soil biogeochemical cycles.

Scatter-hoarding birds disperse seeds to sites unfavorable for plant regeneration.

Scatter-hoarding birds provide effective long-distance seed dispersal for plants. Transporting seeds far promotes population spread, colonization of new areas, and connectivity between populations. However, whether seeds transported over long distances are deposited in habitats favorable to plant regeneration has rarely been investigated, mainly due to methodological constraints. To investigate dispersal patterns and distances of Swiss stone pine (Pinus cembra) seeds we utilized advances in tracking technology to track the movements of their sole disperser, the spotted nutcracker (Nucifraga caryocatactes). We found routine individual movements between single seed harvesting and seed caching site. Harvesting sites of individual birds overlapped, whereas seed caching sites were separated and located on average 5.3 km away from the harvesting site. Interestingly, most distant caching sites were located at low elevations and in spruce forest, where Swiss stone pine does not naturally occur. This suggests that nutcrackers disperse seeds over long distances but that a large portion of these seeds are cached outside the known pine habitat. Therefore, we conclude that the implications of such long-distance seed dispersal movements for plant populations should be carefully considered in combination with the effects of habitat quality on plant recruitment.

Identification and expression of functionally conserved circadian clock genes in lichen-forming fungi.

Lichen-forming fungi establish stable symbioses with green algae or cyanobacteria. Many species have broad distributions, both in geographic and ecological space, making them ideal subjects to study organism-environment interactions. However, little is known about the specific mechanisms that contribute to environmental adaptation in lichen-forming fungi. The circadian clock provides a well-described mechanism that contributes to regional adaptation across a variety of species, including fungi. Here, we identify the putative circadian clock components in phylogenetically divergent lichen-forming fungi. The core circadian genes (frq, wc-1, wc-2, frh) are present across the Fungi, including 31 lichen-forming species, and their evolutionary trajectories mirror overall fungal evolution. Comparative analyses of the clock genes indicate conserved domain architecture among lichen- and non-lichen-forming taxa. We used RT-qPCR to examine the core circadian loop of two unrelated lichen-forming fungi, Umbilicaria pustulata (Lecanoromycetes) and Dermatocarpon miniatum (Eurotiomycetes), to determine that the putative frq gene is activated in a light-dependent manner similar to the model fungus Neurospora crassa. Together, these results demonstrate that lichen-forming fungi retain functional light-responsive mechanisms, including a functioning circadian clock. Our findings provide a stepping stone into investigating the circadian clock in the lichen symbiosis, e.g. its role in adaptation, and in synchronizing the symbiotic interaction.

A Candidate Gene Cluster for the Bioactive Natural Product Gyrophoric Acid in Lichen-Forming Fungi.

Natural products of lichen-forming fungi are structurally diverse and have a variety of medicinal properties. Despite this, they have limited implementation in industry mostly because the corresponding genes are unknown for most of their natural products. Here, we implement a long-read sequencing and bioinformatic approach to identify the putative biosynthetic gene cluster of the bioactive natural product gyrophoric acid (GA). Using 15 high-quality genomes representing nine GA-producing species of the lichen-forming fungal genus Umbilicaria, we identify the most likely GA cluster and investigate the cluster gene organization and composition across the nine species. Our results show that GA clusters are promiscuous within Umbilicaria, and only three genes are conserved across species, including the polyketide synthase (PKS) gene. In addition, our results suggest that the same cluster codes for different, but structurally similar compounds, namely, GA, umbilicaric-, and hiascic acid, bringing new evidence that lichen metabolite diversity is also generated through regulatory mechanisms at the molecular level. Ours is the first study to identify the most likely GA cluster and, thus, provides essential information to open new avenues for biotechnological approaches to producing and modifying GA and similar lichen-derived compounds. GA PKS is the first tridepside PKS to be identified. IMPORTANCE The implementation of natural products in the pharmaceutical industry relies on the possibility of modifying the natural product (NP) pathway to optimize yields and pharmacological effects. Characterization of genes and pathways underlying natural product biosynthesis is a major bottleneck for exploiting the medicinal properties of the natural products. Genome mining is a promising and relatively cost- and time-effective approach to utilize unexplored NP resources for drug discovery. In this study, we identify the most likely gene cluster for the lichen-forming fungal depside gyrophoric acid in nine Umbilicaria species. This compound shows cytotoxic and antiproliferative properties against several cancer cell lines and is also a broad-spectrum antimicrobial agent. This information paves the way for generating GA analogs with modified properties by selective activation/deactivation of genes.

Virus diversity in metagenomes of a lichen symbiosis (Umbilicaria phaea): complete viral genomes, putative hosts and elevational distributions.

Viruses can play critical roles in symbioses by initiating horizontal gene transfer, affecting host phenotypes, or expanding their host's ecological niche. However, knowledge of viral diversity and distribution in symbiotic organisms remains elusive. Here we use deep-sequenced metagenomic DNA (PacBio Sequel II; two individuals), paired with a population genomics approach (Pool-seq; 11 populations, 550 individuals) to understand viral distributions in the lichen Umbilicaria phaea. We assess (i) viral diversity in lichen thalli, (ii) putative viral hosts (fungi, algae, bacteria) and (iii) viral distributions along two replicated elevation gradients. We identified five novel viruses, showing 28%-40% amino acid identity to known viruses. They tentatively belong to the families Caulimoviridae, Myoviridae, Podoviridae and Siphoviridae. Our analysis suggests that the Caulimovirus is associated with green algal photobionts (Trebouxia) of the lichen, and the remaining viruses with bacterial hosts. We did not detect viral sequences in the mycobiont. Caulimovirus abundance decreased with increasing elevation, a pattern reflected by a specific algal lineage hosting this virus. Bacteriophages showed population-specific patterns. Our work provides the first comprehensive insights into viruses associated with a lichen holobiont and suggests an interplay of viral hosts and environment in structuring viral distributions.

Environmental context determines the limiting demographic processes for plant recruitment across a species' elevational range.

Plant recruitment is a multi-stage process determining population dynamics and species distributions. Still, we have limited understanding of how the successive demographic processes depend on the environmental context across species' distributional ranges. We conducted a large-scale transplant experiment to study recruitment of Pinus cembra over six years. We quantified the effects of environmental conditions on four demographic processes and identified the most limiting across and beyond the pines' elevational range over several years. Realized transition probabilities of the demographic processes varied substantially across the species' distributional range. Seed deposition decreased from the lower to the upper elevational range margin by 90%, but this reduction was offset by increased seed germination and seedling survival. Dispersal limitation at the upper range margin potentially stems from unsuitable seed caching conditions for the animal seed disperser, whereas increased seed germination might result from enemy escape from fungal pathogens and favourable abiotic conditions at the upper range margin. Our multi-year experiment demonstrates that environmental context is decisive for the local relevance of particular demographic processes. We conclude that experimental studies identifying the limiting demographic processes controlling species distributions are key for projecting future range dynamics of plants.

The origin of the serpentine endemic Minuartia laricifolia subsp. ophiolitica by vicariance and competitive exclusion.

Serpentine soils harbour a unique flora that is rich in endemics. We examined the evolution of serpentine endemism in Minuartia laricifolia, which has two ecologically distinct subspecies with disjunct distributions: subsp. laricifolia on siliceous rocks in the western Alps and eastern Pyrenees and subsp. ophiolitica on serpentine in the northern Apennines. We analysed AFLPs and chloroplast sequences from 30 populations to examine their relationships and how their current distributions and ecologies were influenced by Quaternary climatic changes. Minuartia laricifolia was divided into four groups with a BAPS cluster analysis of the AFLP data, one group consisted only of subsp. ophiolitica, while three groups were found within subsp. laricifolia: Maritime Alps, north-western Alps and central Alps. The same groups were recovered in a neighbour-joining tree, although subsp. ophiolitica was nested within the Maritime Alps group of subsp. laricifolia. Subspecies ophiolitica contained three different chloroplast haplotypes, which were also found in the Maritime Alps group of subsp. laricifolia. Given its high genetic diversity, subsp. ophiolitica appears to have arisen from subsp. laricifolia by vicariance instead of by long-distance dispersal. Genetic and geographic evidence point to the Maritime Alps populations of subsp. laricifolia as the closest relatives of subsp. ophiolitica. We hypothesize that M. laricifolia was also able to grow on nonserpentine rocks in the northern Apennines during glacial periods when the vegetation was more open, but that only the serpentine-adapted populations were able to persist until the present due to their competitive exclusion from more favourable habitats.